Orbital order, stacking defects and spin-fluctuations in the p-electron molecular solid 
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We examine magnon and orbiton behavior in localized O2 anti-bonding molecular 7r* orbitals using 
an effective Kugel-Khomskii Hamiltonian derived from a two band Hubbard model with hopping 
parameters taken from ab initio density functional calculations. The considerable difference between 
intraband and interband hoppings leads to a strong coupling between the spin wave dispersion 
and the orbital ground state, providing a straightforward way of experimentally determining the 
orbital ground state from the measured magnon dispersion. The near degeneracy of different orbital 
ordered states leads to stacking defects which further modulate spin-fluctuation spectra. Proliferaion 
of orbital domains disrupts long-range magnetic order, thus causing a significant reduction in the 
observed Neel temperature. 

PACS numbers: 71.10.Fd, 75.10.Jm, 75.25.Dk, 75.50.Xx 



Correlated systems have generated considerable inter- 
est in the literature in recent years. The discovery of 
high temperature superconductors and the subsequent 
development and application of correlated methods like 
LDA+U and DMFT 0], has led to remarkable success 
in dealing with strongly correlated systems. At integer 
filling, strongly correlated systems are typically insulat- 
ing and often show antiferromagnetic behavior arising 
from exchange or superexchange processes. Such systems 
include the undoped cuprates, where there is one hole per 
site that can hop in a square lattice of Cu d x 2_ y 2 orbitals, 
and heavy fermion materials like CeCuIns where the Ce 
4/ orbitals weakly couple to the valence states. Multi- 
band correlated systems can also show orbital ordering; 
the earliest successful application of LDA+U found or- 
bital ordering in the KCuF 3 system Q. 

Since correlated behavior is typically the domain of 
materials with 3d and 4/ orbitals, comparatively little 
attention has been given to the study of correlated be- 
havior in p orbital systems. However, local moment mag- 
netism in 2p orbitals has been implicated in several sys- 
tems, such as at polar oxide vacancies [3] and substitu- 
tionals 5]. The occurrence of 2p orbital moments at p- 
type LaAlOa/SrTiOa interfaces [6| is still the only viable 
explanation of the insulating character that is observed in 
these interfaces, where the electron count would suggest 
metallic interface states. Alkali hyperoxides, to be dis- 
cussed below, comprise another likely example. Recent 
calculations [7| suggest that doping of d° (no d electrons) 
magnetic systems can stabilize or even enhance 2p mag- 
netic moments in systems such as ZnO nanowires Q. 

Recently, there has been significant interest in studying 
correlations in solid molecular systems, such as SrN 
which consists of Sr octahedra containing either isolated 
N atoms or N2 dimers, with calculations predicting that 
the magnetic moment is strongly confined to the anionic 
dimers. Calculations on the Rb40 6 system Q and 
Cs 4 6 suggest that these systems would be a half- 



metallic ferromagnets particularly useful for spintronic 
applications, due to the reduced spin-orbit interaction in 
p orbitals. However, these calculations were done within 
weakly correlated density functional theory; more recent 
calculations using LDA+U ll| suggest that the valence 
charge separates to give a mixture of magnetic hyper- 
oxide anions and nonmagnetic peroxide 0\~ anions, 
and an insulating ground state. There seems to be some 
experimental disagreement as to whether Rb406 is con- 
ducting or insulating 11 1. In Rb40g the three dif- 
ferent orientations of the O2 dimers along the principle 
axes of the crystal give rise to frustration of the magnetic 
order. 

Our interest here is in the alkali h ype roxides, taking 
RbC>2 as a specific example. Solovyev [13J has provided a 
study of the sister compound KO2, considering the large 
spin-orbit coupling (SOC) limit. We present here an al- 
ternative viewpoint for Rb02, based on the supposition 
that SOC is not so large an effect, so orbital moments are 
quenched by the crystal field, making the conventional 
redip x ,p v orbitals the natural basis for studying spin and 
orbital phenomena in Rb02- Experimental measurement 
of the Lande g-factor 14[ yields values close to 2, indi- 
cating a mostly spin moment, rather than the g — 4/3 
value expected for large spin orbit coupling. 

The MO2 systems (M=Li,Na,Rb,Cs) exhibit complex 
phase diagrams and low temperature antiferromagnetism 
|15| . The phase diagrams at low temperature consist 
of several structural changes which are minor symmetry 
lowering distortions from the room temperature (aver- 
aged) tetragonal phase which is comparable to a distorted 
rock salt structure with O2" ions playing the role of the 
anion, with the molecular axis pointing along the c direc- 
tion, shown in Fig. [TJd- The Jahn- Teller effect causes the 
O2 molecules to tilt away from the tetragonal axis, an 
effect which is difficult to reproduce in a non-magnetic 
LDA calculation. Below 194 K, Rb02 shows incommen- 
surate superstructure and a mixture of pseudotetragonal, 
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FIG. 1. (Color online) (a) Paramagnetic band plot of Rb02 
showing bands at the Fermi level from both DFT(LDA) and 
tight binding. The bands are filled such that there is one 
hole per site, (b) Conventional unit cell of the tetragonal 
phase of Rb02- (c) Schematic showing tight binding hopping 
parameters in meV. Circles represent OJ anions in a plane 
perpendicular to the molecular axis, triangles represent O^T 
anions in the nearest planes above or below, Rb not shown, 
(d) Isosurface of the tt* Wannier function in its local environ- 
ment, (b) and (d) produced with Vesta 



lows tight-binding hopping parameters to be calculated 
directly. We consider only the four most relevant hop- 
pings, two for nearest neighbors (t xx = (0x\H\H.ix) and 
t X y = (0x\H\Riy), where Ri = i(a,a,c)) and two for 
second neighbors in the plane (t a = (0x\H\R2 X x) and 
tit = (0x\H[R2 y x) with R,2x = ax and R2 y = ay). The 
LDA spin-unpolarizcd band structure and our tight bind- 
ing band structure with four parameters are shown in 
Fig. QJi. A schematic for the hopping channels used in 
the tight-binding model is shown in Fig. [TJ; with the 
numerical values of the hoppings. Note that the largest 
hopping is t xy , which represents nearest neighbor hop- 
ping from \x) to \y) orbitals. Second neighbor hoppings 
from \x) to \y) are forbidden by symmetry. There are 8 
first neighbors with t xy hopping, compared to 2 each for 
t a and t n , and since 8\t xy \ » 2\t a \ + 2|i 7r |, this would 
suggest that the strongest AFM coupling is between near- 
est neighbors. Previous neutron diffraction studies [2(| 
on KO2 show magnetic ordering consistent with nearest 
neighbor spin ordering. 

The non- interacting Hamiltonian is 2x2, with values 

H xx = -8txxl(k) - 2t a cos k x - 2t 7T cos k y (1) 
H yy = -8*ax7(k) - 2U cos k x - 2t a cos k y (2) 
H % = -4t xy cos (|fc z ) sin (^k x ) sin (^k y ) (3) 

where we define the nearest neighbor structure factor 



orthorhombic and monoclinic crystal structures that are 
all slight distortions of the tetragonal phase [IH . 

The standard LDA calculation for RbC>2 produces a 
half-metallic (ferromagnetic) state in the averaged unit 
cell 13|, or an antiferromagnetic metal in a Rb204 super- 
cell, in contrast to experimental reports that MO2 com- 
pounds are insulating [12J. As seen in Fig. QJi, the O2 tt* 
bands near the Fermi level have a bandwidth of about 
1 eV and are well separated from other bands. These 
bands contain 3 electrons per 0.7 ion, so the occupations 



of the tt* and tt* (hereafter 



tt* and \y) = tt*) or- 



bitals are frustrated and likely related to the structural 
transitions, as well as the Jahn- Teller distortion that tilts 
the ions. Since the relevant bands are so narrow and 
the system is Mott insulating, one might consider the 
use of the LDA+U method but since the typical im- 
plementation of LDA+U in DFT codes uses interactions 
between onsite atomic orbitals and the appropriate in- 
teractions in MO2 would be between molecular orbitals, 
a straightforward application of LDA+U fails to produce 
the correct ground state 

to construct a tight- 
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We used the FPLO code 
binding Hamiltonian by projecting Wannier functions 
using symmetry projected orbitals [lil Il9l | correspond- 
ing to the two O2 molecular tt* orbitals in the prim- 
itve cell in a paramagnetic LDA calculation. This al- 



7(k) = cos {\k x ) cos {\k y ) cos (|&z) • (4) 
The interacting Hamiltonian is 

~ (5) 



ia,jf3 



where i,j run over sites and a, (3 € [x,y]. We ap- 
ply second-order perturbation theory to H to get an 
effective Kugel-Khomskii (KK) [2l[ type Hamiltonian 
Hkk — H K n K + Hjzk for nearest and next nearest neigh- 
bor interactions, 
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J a = 4i£/E7, = 4i%/U, J s — Jct + J-k where terms with 
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FIG. 2. (Color online) (a) Spin wave spectrum for Rb02 for 
p-type ordering for both stackings as described in text (black 
solid, ABAB stacking and red dashed, ABCD stacking), and 
spin wave spectrum for xy ordering (blue dotted). Also, the 
nearly dispersionless orbital spectrum is shown for p-type or- 
derings (black dash-dot line), (b) Spin wave spectrum for 
ferroorbital ordering, which is orthorhombic symmetry. 



a single r operators have been neglected. The parame- 
ter Jd = ±( J a — J ff ) depends on the direction of the bond 
between i and j. The r operators can be represented as 
Pauli matrices that operate in orbital space. For numer- 
ical results, we select U = 3 eV, which is equivalent to 
U c s — U — J for intraorbital U and Hund's exchange J 
found previously [l3| for the O2 it* orbitals in KO2 via 
the constrained LDA method This gives J xx = 1.54 
meV, J xy = 4.8 meV, 3 a = 3.89 meV, and J w = 0.34 
meV. 

Unlike previous studies jUH^j on the KK Hamilto- 
nian, which contained only a single parameter t for both 
intraorbital and interorbital hopping, we consider here 
a KK-type Hamiltonian where the hoppings between or- 
bitals channels arc significantly different than hoppings 
within an channel, as shown in Fig. [TJ One previ- 
ous study [26] concluded that due to unusual symme- 
tries present, the KK Hamiltonian could not describe the 
observed order and gapped excitations. In the present 
work, the broken symmetry of the hoppings in this model 
should avoid these difficulties. These hoppings will lead 
to the orbital ground state having a significant impact 
on the spin wave dispersion, which should be measurable 
in experiment. 

To proceed with spin/orbital wave theory a reference 
ground state for both the spin and orbital systems must 
be chosen. From here on, we will consider the electronic 
structure from the hole perspective, so that there is a 



single hole per site. We restrict our magnetic order to be 
antiferromagnetic between with nearest neighbors. 

The ground state orbital ordering for this model is 
where antiferroorbital (AFO) ordering occurs in planes 
so that a given site with, say, \x) occupied would have 
second neighbors (first neighbors in the plane) with \y) 
occupied. We refer to this ordering as p-type. This or- 
dering frustrates the first neighbor orbitals, so alternate 
stackings of the planes will be very close in energy and 
may be degenerate. For these orderings, the spin wave 
dispersion w£ and orbital wave dispersion are given 
by 



(8J) 2 



V^-(^/72(k)) 



[4J XJ/7ro (k) + 4J ra7 „(k)] 2 (8) 

(9) 



where Ji = V J a J-rr and 72 (k) = i(cos fc x +cos k y ). where 
the effect of different stackings is contained in the struc- 
ture factors j m (k) and 7«(k). First we consider the case 
of ABAB stacking, where any given site has the same 
orbital occupation as the sites at displacements (0, 0, c) 
and (0, 0, — c) from it. This results in structure factors 
given by 



7m (k) ABAB = cos {\(k x + k y )) cos(ifc z ) 
7 „(k) ABAB = cos {\{k x - k y )) cos(ifc 2 ). 



(10) 



An alternate stacking, ABCD, where each site has the 
opposite orbital occupation as the sites above and below 
in the z direction, has structure factors of 

7m (k) ABCD = cos % cos % COS % ' 
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7n(k 



\ ABCD 



i sin sin -j- sin 4f 
cos cos cos 4f - 
isin^sin^-sin^. 



(11) 



The dispersions are depicted in Fig. [5^. 

An alternate orbital ordering we consider is where the 
ordering within planes are ferroorbital (FO) ordering. 
This ordering is not stable to orbital fluctuations within 
our KK model, however it may be stabilized in Rb02 by 
effects not considered here, such as the Jahn- Teller dis- 
tortion that causes the canting of the O2 molecules. This 
ordering is of particular interest because if the stacking of 
planes is AFO, then the nearest neighbor spin exchange 
is maximized. With this orbital ordering (hereafter re- 
ferred to as xy ordering), the spin sublattices have differ- 
ent dispersions due to the swapping of J a and J n in each 
plane. The Kugel-Khomskii Hamiltonian yields a spin 
wave dispersion of 



Al-(8J xyl (k)y±A d 



(12) 



with A rn = 8J xy 

(Ja - Jtt) [C0s(fc x ) 



J s [2 — cos(fc x ) — cos(fcy)] and Ad = 
- cos(ky)], and is shown in Fig. [2^. 
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Thus far, the orderings considered all preserve tetrago- 
nal symmetry, but orthorhombic and monoclinic low tem- 
perature phases of Rb02 exist. The final orbital ordering 
we consider is FO ordering for every site in the crystal, 
hereafter called xx ordering, which is orthorhombic sym 
metry. This has spin-wave dispersion ui%. x = yvf 
where — 8J XX + \{J S ~ J<? cos k x — J n cosk y ) and 
£>k = 8J xx 7(k). Again, this configuration is not stable 
with respect to orbital fluctuations, but it may be sta- 
bilized by effects not considered here, which a study of 
KO2 indicates is the case 27J. The magnon dispersions 
is shown in Fig. [2)3. 

Orbital Order 
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p-ABAB 


p-ABCD 
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-12.680 


-12.680 


-14.795 


-14.795 


E s 


-0.411 


-3.266 


-2.333 


-2.288 


Eo 








-0.010 


-0.010 


Etot 


-13.091 


-15.946 


-17.138 


-17.093 



TABLE I. Energies contributing to the ground state energy 
in meV. -Bo is the classical ground state energy, Es{Eo) is 
the quantum correction from spin(orbital)-wave theory. 

The ground state energies for these four orbital or- 
derings arc listed in table [H where we find that the 
ABAB stacking of planar orbital ordering is the lowest 
energy. For an average nearest neighbor spin exchange 
of J = 3.17 meV (in p-type orderings), high temperature 
series expansion [28| would predict T N f» 1.4 J = 51 K, 
much higher than the observed T^r = 15 K. The frus- 
trated exchanges within planes would reduce this value 
somewhat, but not enough to give a prediction reason- 
ably close to the experimental transition. 

As expected, the two stackings examined for the pla- 
nar orbital ordering are very nearly degenerate, with only 
quantum fluctuations in the spin waves breaking the de- 
generacy at this level of approximation. It is rather 
clear that due to the strong asymmetry between hopping 
within an orbital channel and hopping between orbital 
channels, the orbital ground state significantly impacts 
the spin wave spectrum and could be inferred from a mea- 
surement of the low temperature magnons. The planar 
orbital orderings don't significantly impact the spin wave 
dispersion, so even if the stacking is disordered magnon 
excitations should be coherent. Energetically the next 
state above p- type ordering is the xy ordering, which 
is ^1.2 meV « 14 K higher in energy than p-type or- 
dering. Above this temperature orbital domains should 
proliferate. The strong modulation of exchange constants 
at the domain boundaries should nucleate magnetic do- 
mains , leading to the low observed Neel temperature 
of Rb0 2 of 15 K. 

Orbitons are quite difficult to measure experimentally. 
They do not couple directly to neutrons, the standard 



measurement technique for magnons. Recently orbitons 
have been measured in titanates via X-rays [3fJ , but the 
inference of the orbital dispersion is very indirect. The 
measured dispersion is negligible, suggesting [3(| that X- 
rays modulate bonds resulting in a much bigger scatter- 
ing from two-orbitons than from single orbitons. It may 
be that in narrow band molecular oxide systems the X- 
rays will create single orbiton excitations with all bond 
modulation inside a unit cell. If that is the case, the 
X-ray measured dispersion in RbC>2 should be clearly 
observable, if a similar experiment can be done. Ne- 
glected here is the Jahn- Teller effect which rotates the 
O2 molecules and breaks the orbital degeneracy, which 
will select a particular ground state orbital ordering. 

We have examined independent magnon and orbiton 
excitations in RbC>2 within spin/orbital wave theory, 
finding considerable coupling between the easily mea- 
sured magnon excitations and the difficult to measure 
orbital ground state. This strong coupling arises from 
the large anisotropy in the hopping parameters in the 
O2 7r* bands. This indicates that MO2 materials are 
attractive for studying the interplay between orbital or- 
dering and p-electron magnetism. E. R. Y. and W. E. P 
were supported by DOE SciDAC Grant No. DE-FC02- 
06ER25794. The authors would like to thank J. Kunes, 
R. T. Scalettar, and C. Felser for stimulating conversa- 
tion. 
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